A-baryon electromagnetic form factors in lattice QCD 
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We develop techniques to calculate the four A electromagnetic form factors using lattice QCD, 
with particular emphasis on the sub-dominant electric quadrupole form factor that probes defor- 
mation of the A. Results are presented for pion masses down to approximately 350 MeV for three 
cases: quenched QCD, two flavors of dynamical Wilson quarks, and three flavors of quarks described 
by a mixed action combining domain wall valence quarks and dynamical staggered sea quarks. The 
magnetic moment of the A is chirally extrapolated to the physical point and the A charge density 
distributions are discussed. 
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Lattice Quantum Chromodynamics (QCD) provides a 
well-defined framework to directly calculate hadron form 
factors from the fundamental theory of strong interac- 
tions. Form factors characterize the internal structure 
of hadrons, including their magnetic moment, their size, 
and their charge density distribution. Since the A(1232) 
decays strongly, experiments [l], 0] to measure its form 
factors are harder and yield less precise results than for 
nucleons [3j, |J] . In this work, we compute A form factors 
using lattice QCD more accurately than can be currently 
obtained from experiment. 

A primary motivation for this work is to understand 
the role of deformation in baryon structure: whether any 
of the low-lying baryons have deformed intrinsic states 
and if so, why. Thus, a major achievement of this 
work is the development of lattice methods with suffi- 
cient precision to show, for the first time, that the elec- 
tric quadrupole form factor is non-zero and hence the A 
has a non-vanishing quadrupole moment and an associ- 
ated deformed shape. Unlike the A, the spin-1/2 nucleon 
cannot have a quadrupole moment, so the experiment of 
choice to explore its deformation has been measurement 
of the nucleon to A electric and Coulomb quadrupole 
transition form factors. Major experiments 0] have 
shown that these transition form factors are indeed non- 
zero, confirming the presence of deformation in either the 
nucleon, A, or both [H, and lattice QCD yields com- 
parable non-zero results [loL [Tl| . Our new calculation of 
the A quadrupole form factor, coupled with the nucleon 
to A transition form factors, should in turn shed light on 
the deformation of the nucleon. 

In order to evaluate the A electromagnetic (EM) form 
factors to the required accuracy, we isolate the two 
dominant form factors and the sub-dominant electric 
quadrupole form factor. This is particularly crucial for 
the latter since it can be extracted with greater precision, 
although it increases the computational cost. Our tech- 
niques are first tested in quenched QCD We then 
calculate form factors using two degenerate flavors of dy- 



namical Wilson fermions, denoted by Np — 2, with pion 
masses in the range of 700 MeV to 380 MeV 0J] . Fi- 
nally, we use a mixed action with chirally symmetric do- 
main wall valence quarks and staggered sea quarks with 
two degenerate light flavors and one strange flavor [l5j . 
denoted by Np = 2+1, at a pion mass of 353 MeV. 
Using the results obtained with dynamical quarks, we 
extrapolate the magnetic moment to the physical point. 
We extract the quark charge distributions in the A, and 
discuss their quadrupole moment. 

The A matrix element (A(p f ,s f )\j^ M \A(p i ,s i )), 
where jg M is the electromagnetic current, can be 
parametrized in terms of four multipolc form factors that 
depend only on the momentum transfer q 2 = —Q 2 = 
(pf —pi) 2 [la]. The decomposition for the on shell 7*AA 
matrix element is given by 

(A( Pf ,s f )\j£ M \A(p i ,s i )) = A u«r(p/,s/)O a/1T u r (p i ,s i ) 
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where a\(q 2 ), d2(q 2 ), c\{q 2 ), and C2{q 2 ) are known linear 
combinations of the electric charge form factor GEo(q 2 ), 
the magnetic dipole form factor Gmi(? 2 ), the electric 
quadrupole form factor Gs^o 2 ), and the magnetic oc- 
tupole form factor Gm3(<? 2 ) [1Z|, an d A is a known factor 
depending on the normalization of hadron states. These 
form factors can be extracted from correlation functions 
calculated in lattice QCD [TtJ • We calculate in Euclidean 
time the two- and three-point correlation functions in a 
frame where the final state A is at rest: 
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where j M is the electromagnetic current on the lattice, J 
and J are the A + interpolating fields constructed from 
smeared quarks 0, T 4 = |(1 + 7 4 ), and T fc = iT 4 -f 5 j k . 
The form factors can then be extracted from ratios of 
three- and two-point functions in which unknown nor- 
malization constants and the leading time dependence 
cancel 
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For sufficiently large t/ — t and t — tj, this ratio exhibits a 
plateau R(T, t,q) — > n(r, tf), from which the form factors 
are extracted, and we use the particular combinations 
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The connected part of each combination of three-point 
functions can be calculated efficiently using the method 
of sequential inversions jl8j ]. At present, it is not yet 
computationally feasible to calculate the small correc- 
tions arising from disconnected diagrams. The known 
kinematical coefficients K\, K3, K4 are functions of 
the A mass and energy as well as of \i and q. The combi- 
nations above are chosen such that all possible directions 
of \x and q contribute symmetrically to the form factors at 
a given Q 2 [19]. The over-constrained system of Eqs. ([4j 
[6|) is solved by a least-squares analysis, and Ge2{Q 2 ) can 
also be isolated separately from Eq. ©. 

The details of the simulations are summarized in Ta- 
ble HI In each case, the separation between the final 
and initial time is tt — tj > 1 fm and Gaussian smear- 
ing is applied to both source and sink to produce ade- 
quate plateaus by suppressing contamination from higher 
states having the quantum numbers of the A(1232). For 
the mixed-action calculation, the domain-wall valence 
quark mass was chosen to reproduce the lightest pion 
mass obtained using Np = 2 + 1 improved staggered 
quarks pjiHtJ. 

The results for Geo(Q 2 ) are shown in Fig. Q] as a func- 
tion of Q 2 at the lightest pion mass for each of the three 
actions. For Wilson fermions, we use the conserved lat- 
tice current requiring no renormalization. The local cur- 
rent is used for the mixed action, and the renormaliza- 
tion constant, Zy — 1.0992(32), is determined by the 
condition that Geo(0) equals the charge of the A in 
units of e. As can be seen, all three calculations yield 
consistent results. The momentum dependence of the 
charge form factor is described well by a dipole form 

Geo(Q 2 ) = 1/ U + -jprA ■ To compare the slopes at 



TABLE I: Lattice parameters and results. N con { denotes the 
number of lattice configurations, -J (r 2 ) gives the charge ra- 
dius, is the A + magnetic moment in nuclear magnetons 
and is the A + quadrupole moment. 



N con { 


rrw [GeV] 


mA [GeV] VP) [fm] 


ma+ [m] 


Qf 
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Quenched Wilson, 32 a x 64, a 


= 0.092 fm 


200 
200 
200 


0.563(4) 
0.490(4) 
0.411(4) 


1.470(15) 0.6147(66) 
1.425(16) 0.6329(76) 
1.382(19) 0.6516(87) 


1.720(42) 
1.763(51) 
1.811(69) 


0.96(12) 
0.91(15) 
0.83(21) 


N F = 


2 Wilson, 24 d x 40(32 for lightest pion), a = 


0.077 fm 


185 
157 
200 


0.691(8) 
0.509(8) 
0.384(8) 


1.687(15) 0.5279(61) 
1.559(19) 0.594(10) 
1.395(18) 0.611(17) 


1.462(45) 
1.642(81) 
1.58(11) 


0.80(21) 
0.41(45) 
0.46(35) 


N F 


= 2 + 1, Mixed action, 28^ x 64, 


a = 0.124 fm [21] 


300 


0.353(2) 


1.533(27) 0.641(22) 


1.91(16) 


0.74(68) 
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FIG. 1: The electric charge form factor versus Q 2 . The green 
(red) line and error band show a dipole fit to the mixed action 
(quenched ) results. 

Q 2 = 0, we follow convention and show in table U the so- 
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called "rms radius" [TlJ (r 2 ) 

The momentum dependence of Gmi(Q 2 ) is displayed 
in Fig. O To extract the magnetic moment, an extrapo- 
lation to zero momentum transfer is necessary. Both an 
exponential form, Gmis~® / Ajji , and a dipole describe 
the Q 2 -dependence well, and we adopt the exponential 
form because of its faster decay at large Q 2 , in accord 
with perturbative arguments. The larger spatial volume 
for the quenched and mixed action cases yields smaller 
and more densely spaced values of the lattice momenta 
and correspondingly more precise determination of the 
form factor than for the smaller volume used with dy- 
namical Wilson fermions. In Fig. [21 we show the best 
exponential fit and error band for the mixed action and 
quenched results. As can be seen, results in the quenched 
theory and for Np = 2 Wilson fermions are within the 
error band. The magnetic moment in natural units is 
given by jtiA = Gi,/i(0)e/(2mA), where mA is the A 
mass measured on the lattice and Gmi(0) is from the 
exponential fits. In Table U we give the values of the 
A + magnetic moment in nuclear magnetons e/(2Afjv), 
with Mn the physical nucleon mass. The magnetic 
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FIG. 2: The magnetic dipole form factor. The green (red) line 
and error band show an exponential fit to the mixed action 
(quenched )results. 
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FIG. 3: The magnetic dipole moment in nuclear magnetons. 
The value at the physical pion mass (filled square) is shown 
with statistical and systematic errors The solid and 

dashed curves show the chiral extrapolation and theoretical 
error estimate [2^] . 

moments of the A + and A ++ are accessible to experi- 
ments 0, , which presently suffer from large uncertain- 
ties. The magnetic moment as a function of m% is shown 
in Fig. [31 together with a chiral extrapolation to the phys- 
ical point [22], which lies within the broad error band 
Ma+ = 2.T±\%{stat.) ± 1.5(syst.) ± 3.0(theory)fi N 
The A moments using an approach similar to ours are cal- 
culated only in the quenched approximation [U [H, [HJ • 
Our magnetic moment results agree with recent back- 
ground field calculations using dynamical improved Wil- 
son fermions [2u| . which supersede previous quenched 
background field results [26| , The spatial length L s of 
our lattices satisfies Lgin^ > 4 in all cases except at the 
lightest pion mass with Np — 2 Wilson fermions, for 
which LgHn^ = 3.6. For that point, the magnetic mo- 
ment falls slightly below the error band, consistent with 
the fact that Ref. 25] shows that finite volume effects 
decrease the magnetic moment. 

The electric quadrupole form factor is particularly in- 
teresting because it can be related to the shape of a 
hadron, and lattice calculations for each of the three 
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FIG. 4: The electric quadrupole form factor. The notation is 
the same as that in Fig. [T] 

actions are shown in Fig. 0] with exponential fits for 
the quenched and mixed action cases. Just as the 
electric form factor for a spin 1/2 nuclcon can be ex- 
pressed precisely as the transverse Fourier transform of 
the transverse quar k charge density in the infinite mo- 
mentum frame [27], a proper field-theoretic interpreta- 
tion of the shape of the A(1232) can be obtained by 
considerin g th e quark transverse charge densities in this 
frame [2cl |29L [30 ]. With respect to the direction of the 
average baryon momentum P, the transverse charge den- 
sity in a spin-3/2 state with transverse polarization s± is 
defined as : 
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where the photon transverse momentum q± satisfies q 2 = 
Q 2 , J + = J° + J 3 , and b specifies the quark position in 
the xy-pl&ne relative to the A center of mass. Choos- 
ing the A transverse spin vector along the x-axis, the 
quadrupole moment of this two-dimensional charge dis- 
tribution is defined as fijl: 
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In terms of the A EM form factors [19( , 
Qf = i {2 [Gmi(0) - 3e A ] + [G E2 (0) + 3e A ]} (9) 

The term proportional to [Gmi(0) — 3ca] is an electric 
quadrupole moment induced in the moving frame due 
to the magnetic dipole moment. For a spin-3/2 parti- 
cle without internal structure, Gmi(0) = 3eA, G^2(0) = 
— 3eA [19ll3ll]. and the quadrupole moment of the trans- 
verse charge density vanishes. Hence Q^ ± , and thus the 
deformation of the two dimensional transverse charge 
density, is only sensitive to the anomalous parts of the 
spin-3/2 magnetic dipole and electric quadrupole mo- 
ments, and vanishes for a particle without internal struc- 
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ture. The analogous property holds for a spin-1 parti- 
cle (30| . indicating the generality of this description in 
terms of transverse densities. 




FIG. 5: Quark transverse charge density in a A + polarized 
along the z-axis, with s± — +3/2. The light (dark) regions 
correspond with the largest (smallest) values of the density. 

Fig. [5] shows the transverse density for a A + 

with transverse spin sj_ = +3/2 calculated from the fit 
to the quenched Wilson lattice results for the A form 
factors (which has the smallest statistical errors of the 
three calculations). It is seen that the A + quark charge 
density is elongated along the axis of the spin (prolate) . 
This prolate deformation is robust in the sense that the 
values for Qf obtained from Eq. ((9]) and given in Table U 

t 2 

are all consistently positive. 

In the case of the magnetic octupole form factor [Tgj | . 



which is related to the magnetic octupole moment 0a = 
GW3(0)e/2m A , our statistics are insufficient to distin- 
guish the result from zero. 

In summary, a formalism for the accurate evaluation 
of the A electromagnetic form factors as functions of q 2 
has been developed and used in quenched QCD and full 
QCD with N F = 2 and 2+1 flavors. The charge ra- 
dius and magnetic dipole moment were determined as a 
function of and the dipole moment was chirally ex- 
trapolated to the physical point. The electric quadrupole 
form factor was evaluated for the first time with sufficient 
accuracy to distinguish it from zero. The lattice calcula- 
tions show that the quark density in a A + of transverse 
spin projection +3/2 is elongated along the spin axis. 
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